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Abstract: We examine the dynamics of a self-gravitating magnetized neutron gas as 
a source of a Bianchi I spacetime described by the Kasner metric. The set of Einstein- 
Maxwell field equations can be expressed as a dynamical system in a 4-dimensional phase 
space. Numerical solutions of this system reveal the emergence of a point-like singularity as 
the final evolution state for a large class of physically motivated initial conditions. Besides 
the theoretical interest of studying this source in a fully general relativistic context, the 
resulting idealized model could be helpful in understanding the collapse of local volume 
elements of a neutron gas in the critical conditions that would prevail in the center of a 
compact object. 
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1. Introduction 

Critical stellar configurations, such as white dwarfs, neutron stars, supermassive stars, 
relativistic star clusters and black holes are important astrophysical systems in which rel- 
ativistic effects cannot be ignored. These astrophysical systems, denoted generically by 
the term "compact objects", can be thought of as natural laboratories to understand a 
wealth of phenomena relevant to theoretical and experimental physics under critical condi- 
tions. These conditions are ideally suited to propose and test theoretical models of strong 
magnetic fields associated with self-gravitating sources under critical conditions. In such 
conditions, gravitation couples with other fundamental interactions (strong, weak and elec- 
tromagnetic) and from this interplay important clues of their unification could emerge. 

The presence and effects of strong magnetic fields in compact objects has been studied 
in the literature (see |], ^, ||], |||, [|, ||], and references quoted therein). At a very basic 
level, the simplest approach is to consider various types of self-gravitating and magnetized 
plasmas of neutron matter under equations of state that are appropriate for the critical 
conditions of a compact object [jj], [sj. Since astrophysical systems exhibit, in general, 
angular momentum and pressure anisotropies, it is interesting to examine how these effects 
could influence their stability. 

Following the basic known methodology [[jj , an of equation of state based on a nuclear 
ferromagnetic model (not related to electric currents) was examined in Q, with the purpose 
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of studying the interplay of pressure anisotropy and the magnetic field within a Newtonian 
framework. Hence, in the present work we consider a relativistic generalization that will 
allows us to examine (under specific restriction of the geometry) the evolution of this type 
of magnetized neutron gas under strong gravity. In particular, our aim is to provide the 
minimally basic qualitative and quantitative elements to address the question of whether 
the magnetic field can "slow down" , or even reverse or stop, gravitational collapse, as well 
as the issue of the evolution to a stable configuration. 

It is evident that a proper and comprehensive study of a magnetized fermion g as , as a 
source of a compact object, requires a spacetime with axial symmetry (or without symme- 
tries), leading to dynamical equations that must be solved by means of hydrodynamic codes 
of high complexity g |J]. In previous articles |10| , [llj we studied a self-gravitating 

gas of magnetized electrons, described by an appropriate equation of state, as the source 
of an anisotropic Bianchi-I spacetime described by the Kasner metric. While a Bianchi 
I model is obviously inadequate as the metric of a compact object, it is among the less 
complicated geometries compatible with a magnetized source. Thus, the evolution of such 
a source under a much simplified Bianchi geometry is basically a toy model, but as such 
it can still be useful to obtain qualitative features of the sources under local critical condi- 
tions, specially in conditions near the center of the object where angular momentum plays 
a minor dynamical role. These qualitative results could provide a better understanding, 
and/or a useful approximated description, of the behavior of local internal volume elements 
near the center of a more realistic configuration. In this article we extend previous work 
on the electron gas to the case of a neutron gas. Hence, we follow a similar methodology 
based on re-writing the field and conservation equations for the magnetized neutron gas 
in the Bianchi I geometry as a dynamical system, evolving in a 4-dimensional phase space. 
This system is then analyzed qualitatively and numerically. 

It is worthwhile remarking the basic differences between the magnetized electron gas 
examined in ffl, 1C, 11] and the neutron gas that we consider in this article. Electrons 
interact with a magnetic field through their charge, leading to the Landau diamagnetism 
characterized by a quantization effect associated with the Landau energy levels. However, 
neutrons interact with the magnetic field by their anomalous magnetic moment ( AMM) , in 
the context of Pauli's paramagnetism and the equations of Pauli-Dirac. Consequently, one 
expects the magnetic neutron interaction to be weaker, though conditions of degenerated 
neutron gases in compact objects are also expected to be more critical (because of the higher 
densities) than those of degenerated electron gases. Therefore, relativistic effects of gravity 
are more likely to play a dominant dynamical role in neutron gases. Still, it is important 
to mention that a self-gravitating and magnetized neutron gas is a simplified model of a 
source for a compact object, as protons and chemical equilibrium potentials should also 
be considered to evaluate local effects 0, |§. Nevertheless, the magnetized neutron 
gas already exhibits important qualitative and quantitative differences in comparison with 
electron gases previously examined. 

The paper is organized as follows: we derive in section II the equation of state for 
a magnetized neutron gas. The Einstein-Maxwell field equations for the Kasner metric 
and this source are derived in section III. In section IV we examine the limit of a weak 
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magnetic field. The set of ordinary non-linear differential equations that follows from the 
field equations is transformed into a set of evolution equations in sections V and VI. This 
dynamical system is analyzed qualitatively and numerically in section VII. Our conclusions 
are presented in sections VIII. 



2. Magnetized neutron gas: the equation of state. 

The main properties of gas of magnetized degenerated neutrons are well known (see fl], |9j, 
|j~5||). Considering the grand canonical ensemble, a subsystem can be thought of as a local 
volume of the neutron gas under the influence of a magnetic field H associated with the rest 
of the system (in an astrophysical context this could be good approximation to a suitable 
volume element inside a compact object). Because of this field, the subsystem becomes 
polarized, leading to a magnetization vector that satisfies the relationship: H = B — iirM. 
The field H can be though of as "external" to the subsystem, while B can be "internal" to 
any particle within the subsystem, which feels (in addition to H) the contribution 4irAi 
from particles from the rest of the subsystem. 

The equation of state of the neutron gas follows from calculating the energy spectrum 
of the particles that make up the system. We can obtain this energy spectrum from the 
Dirac equation for neutral particles with anomalous magnetic moment: 

where a^x = ^(7 m 7a — 7a 7^) is the spin tensor, F^ u is the electromagnetic field tensor (we 
have set h = c = 1) and Vl/ is the Dirac field. Solving equation fl2.1|) leads to the following 



energy spectrum |]15[, |lq|, |17|| : 



E n (p, B, ri) = \jp\ + (sjpl +ml + V qB)\ (2.2) 

where p\\,p± are, respectively, the components of momentum in the directions parallel and 
perpendicular to the magnetic field B, m n is the mass of the neutron, q = — 1.91/ijy is 
the neutron magnetic moment {p^ = e/2m p is the nuclear magneton), rj = ±1 are 03 
eigenvalues corresponding to the two possible orientations (parallel or antiparallel) of the 
neutron magnetic moment with respect to the magnetic field. 
The thermodynamical grand potential then takes the form 

n = -kTlnZ, (2.3) 

where k is the Boltzman constant, T is the temperature, Z = Tr(p) is the partition function 
of the system, p = e -( H -^ N )/ kT ; jg the Hamiltonian, p is the chemical potential and 
TV" is the number of particles operator. 

The energy-momentum tensor associated with an external constant magnetic field 
takes the form: 

r " " = + E ^|V 4 S% + 4i^iV^ - n 6%, (2.4) 
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so that in the limit of zero magnetic field limit we obtain the perfect fluid tensor: v = 
pd^ v — (p+ U)5^ A 5 A U . The components of the tensor fl2.4|) are: 

r 3 3 = = v\\ = -n = p, (2.5) 
r\ = r 2 2 =p± = -n-BM = p -bm, (2.6) 

T 4 4 = -U = -TS -fiN-Q, (2.7) 

where S is the specific entropy, N = —dVt/d^i is the particle number density, A4 = —dQ/dB 
is the magnetization, U is the energy density and p is the pressure in the direction of the 
magnetic field. 

The thermodynamical potential can be split in the following two parts: 

where the first term in the right hand side is the statistical contribution and the second is 
the vacuum contribution M. Explicitly, we have 



1 f 00 

"-j-27 E / P±d P± dpMf + (»n,Or(»n,0), (2-8) 



where £ = 1/ksT and / ± (/^ n) C) = (l + e^ nT/in ^) represent, respectively, the contributions 
from the particles and from the antiparticles. The vacuum term is given by the expression: 



ttvn = -r^-z }~ \ P±dp±dp 3 E n , (2.9) 

which is divergent, but can be renormalized, and for fields of intensity B < 10 18 G its 
contribution is not important |J, hence we will neglect this term in the remaining of the 
present article. 

Equation ( |2.8j ) can be easily integrated in the case that concerns us (T = 0), leading 

to 



Q sn — —A 

T7=l, — 1 



Vf* (1 + 7/^(5^ -3)^ (1 + y?/3) 3 (3- 77/?) T V^ 3 

i ~ 1 ~ J^t] 7, — s r) 



12 24 24 

where we have introduced the following expressions: 



(2.10) 



fv = vV 2 - (1 + VP) 2 , s v = 7i~ arcsm ( 1 + 71 ) , A* = — (2.11) 

2 \ fj, J m n 



L v = ln £f4 , /3 = ^, (2.12) 



with B c = m n /q ~ 1.56 x 10 20 G being the critical field and A = m^/47r 2 /i 3 c 3 = 4.11 x 
10 36 erg cm -3 . 

All thermodynamical variables of the system follow readily from the thermodynamical 
potential $7. In particular, by computing neutron density and magnetization in this way 
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we obtain N = NqTn, M. = M.qTm, where iVo = A/m n , Mq = Noq, while the coefficients 
^n,Fm take the form: 



T7=l, — 1 



r//3(l + 7?/?)/, 7,/3/x 2 

■7T H ~ L 



At 



17=1,-1 



{l-2nP)nf n (l + r?/3) 2 (l- 77/3/2) 



^7 g 



Therefore, given (Q) and ( 2.1 0| ) , we can write the equation of state for a relativistic 
degenerated neutron gas in the presence of an external magnetic field as: 



u = fi n N + n = \r u ((3,n), 
p = -n = \r P (p,n), 
m = bm = \pr M (P,n), 



(2.13) 
(2.14) 
(2.15) 



where 



E 



12 24 



(1+7^)3(3-7^) 

24 



/jr^ - r P . 



Let us remark at this point that in ( 2,lC| )-( p,15D we are summing over the magnetic 
moments parallel (rj = —1) or antiparallel (rj = 1) to the magnetic field (i.e. the well 
known Pauli Paramagnetism). The choice rj = dbl, [|| is equivalent to consider different 
phases of the system. The appearance of the threshold of the value of the magnetic field 
for each one of these cases can be seen if we analyze the expressions of the functions 
and s v in ( 2.11]) . We have take /3 > therefore starting with f v : 



however as 77 



> (l + ??/3) 2 = l + 2r ? /3 + r ? 2 /3 2 . 
-1,1 then always rj 1 = 1 and we can rewrite: 

//>(r/ + /3) 2 . 



(2.16) 
(2.17) 



hence; 

03 + 77-/x)(/? + 7 / - + /x)<O. (2.18) 

Now we have two possibilities. But, it is easy to realize that the only acceptable 
possibility is: 

-/x-77 < p <fi-ri. (2.19) 
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This is exactly the restriction that comes from the function s^: 

.1 + 77/3 



< 1. (2.20) 



Therefore, from either one (2.1£) or ( |2.20| ) the following constraints are obtained for the 
magnetic field: 

if v = l^> -1-H<P<H-1, (2.21) 
if t] = -1=> l-n<P<l + n. (2.22) 

These last inequalities restrict the values of the magnetic field. Thus neutrons with AMM 
aligned to the magnetic field have only magnetic field values in the range between 1 — \i < 
P < 1 + (j,. Similarly, for neutrons with AMM oriented antiparallel to the magnetic field 
can take magnetic field values in the interval — 1 — /x</3</x — 1. In the particular case 
fj, = 1 we have. 

if ? ? = -1, < < 2, (2.23) 
if rj = l, -2 < p < 0. (2.24) 

This means that for fields smaller (or equal) than two times the critical magnetic field, 
the neutrons are aligned with r] = — 1 (the AMM is parallel to the magnetic field). Hence, 
neutrons with rj = 1 are forced to invert its sense. 

3. Einstein-Maxwell equations. 

If a local self-gravitating volume of a magnetized neutron gas evolves in the conditions 
prevailing inside a compact object, relativistic effects will be important. This implies 
that local dynamics must be studied in the framework of General Relativity, by means of 
Einstein's field equations: 

together with energy-momentum balance and Maxwell equations: 

T^- v = 0, (3.2) 
F^. v = 0, F [fMU . a] = 0. (3.3) 

where k = 87tGn, with Gn being Newtons's gravitational constant, while square brackets 
denote anti-symmetrization in [iv\ a. The energy-momentum tensor T% associated to 
the magnetized neutron gas is given by (||), with the relevant thermodynamical potentials 
obtained through Statistical Mechanics and complying with the appropriate equation of 
state (as discussed in the previous section). This tensor can also be given in terms of the 
4- velocity u a field as: 

T a p = (U + P)u a u p + P5% + U a p, P = p- (3.4) 
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We will consider the field equations (|3.1|) -(p^) with (|3.4|) as the source tensor of a 
Bianchi I model described in the coordinate representation known as the Kasner metric: 

ds 2 = Qtitfdx 2 + Q 2 (t) 2 dy 2 + Qz{t) 2 dz 2 - dt 2 . (3.5) 

which suggests considering a comoving geodesic 4-velocity u a = 8f, so that the anisotropic 
pressure tensor in (plf ) in the coordinates [x, y, z, t] takes the form: 

n a /3 = dia g [n,n,-2n,o], n = -^, n% = o. (3.6) 

The field equations (3.1) for (|3.4|) and (|3.6|) take the form: 



n x Q2Q3 Q2 Q3 , nKA\ dn\ 

-G x = 7777- + 7T + 7T = ~ k{ $ ~ BM >> 

~ GV y = 7T7T + 7T + 7T = "(p - S-M), (3.8) 

Q1Q2 . Qi . Q2 /„ nX 

- G2 = Q 1 ^ + gr + ^ = ;^' (3 - 9) 

Q1Q2 Q1Q3 Q2Q3 



where A = A- a u a = Aj- Energy-momentum balance ( |3.2| ) leads to: 

U = ^(p + U)-(^ + ^){-BM +p + U). (3.11) 

W3 Vl V2 

while the Maxwell equations equations (|3.3D imply: 



Qi Q2 IB , 



The Einstein-Maxwell equations ([T7|, 3.11 and |3,12] ) are non-linear second order ordinary 



differential equations for the metric functions Qi,Q 2 ,Q3 an d U. In order to treat this 
system numerically, it is necessary to introduce a new set of variables that will transform 
it into a system of first order evolution equations. However, before undertaking this task, 
we examine in the following section the weak field limit. 



4. Limit of weak magnetic field. 

The discussion of the limit of weak magnetic field is important in to illustrate the connection 
between our quantum magnetic field and a classic Maxwellian field in the context of a 
magnetohydrodynamic treatment. In sections 2 we have shown equations strongly linked 
with quantum magnetic field. 

It is important to emphasize that the term "quantum magnetic field" involves the semi- 
classical interaction between the magnetic field and the anomalous magnetic moment. This 
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approach implies a theoretical connection between the equation of state introduced in the 
previous section and a QED framework. As a consistency condition, this framework must 
allow for a classical Maxwellian limit that should arise from a series expansion around the 
zero of the magnetic field (3 = 0. The leading term of this expansion should lead to the 
well known energy- momentum tensor for a Maxwellian magnetic field [21|. In general, this 
type of series expansions can be done as follows: 



k = eO u £r " m - 01/32 + ° {n (41) 

n=0 



oo on 



PH = E(-^T) U ~~J — Pi + as/3 2 + 0(/3 4 ), (4.2) 



n=0 

<x> 



»rr on 

^ = £( aST) U ^ - + ao/3 2 + 0(/3 4 ), (4.3) 

n=0 ^ 

where it is easy to see that: p\ = m{fi = 0) = p±((3 = 0), ao = (d 2 U /d(3 2 ) L =0 /2, ai = 
(d 2 p±/d(3 2 ) \ p=0 /2, a 3 = (9 2 p||/a/3 2 ) |^ =0 /2, C/ = C^(/3 = 0). Hence, all of the former are 
functions only of the dimensionless chemical potential \x. 

As shown in the examples in the literature of classical magnetic fields in the context of 
a Bianchi-I geometry |L9|, |20[] , the leading magnetic field terms are quadratic. This suggest 
to truncate the series in (|4,1|)-( |4~3| ) in the quadratic terms, as the higher order terms, 
like /3 4 , /3 6 , are (in general) extremely small multipole contributions. If we assume 
that the velocity fluctuations of the plasma tend to a zero average macroscopically, and 
that the medium does not undergo any bulk motion, then these contributions can be 
safely neglected (though typically higher velocities could arise from thermal fluctuations or 
quantum disorder). Under these assumptions, the energy-momentum tensor of a neutron 



gas with a minimally coupled magnetic field can always be written in the form [19|: 



Tun, = (U + U mag )u^U v + (po + Pmag)h^ + ■ (4.4) 

where p ma g and U mag are, respectively, the magnetic pressure and the magnetic energy. 

It is important to remark that po in Q4.4Q is the isotropic contribution to the pressure 
of the system, which (in general) can depend on the chemical potential and the magnetic 
field. On the other hand, p\ is the pressure for the case without magnetic field ft = 0. In 
the case under consideration we have: 

T» v = diag[pi - ai/3 2 ,pi - ai (3 2 , Pl + a 3 /3 2 , -U - a /3 2 ]. (4.5) 

Comparison with equations fl4,5D -(fPl) leads to: 

H 2 

Prnag = — = a /3 2 /3, (4.6) 
= ^Y= «0/3 2 - (4.7) 
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where H 2 = H^H^ = 2ao/3 2 and are the components of the magnetic field, which we 
have assumed to point in the z-direction. This assumption is consistent with the fact that a 
small volume element around the center of a compact object is approximately homogeneous 
and the rotation axis furnishes a privileged direction, which we can always align with the 
z axis: 

H a = (0,0,V&?-,0), H a = (0,0,V2^/3Q 3 ,0). (4.8) 

The tensor li^u 9 is then the projected symmetric trace-free tensor representing anisotropic 
pressures that comes from the magnetic field. It can be written as: 

(ir^)t = diag[-i( ai + a 3 )/3 2 , -i( ai + a 3 )/3 2 , |( ai + a 3 )/? 2 , 0]. (4.9) 



Is important to point out that p\ in ( |4.5| ) is, in general, different from po- Thus: 

/3 2 

Po = Pi ~ (ao + 2oi - 03) — - (4.10) 



where pi = P(r=q\ = P[B=o) m ( P-14 ) and (3.4) is the pressure without magnetic field. Only 



when the magnetic field is zero, /3 = 0, the pressures coincide: po = Pi an d the energy- 
momentum tensor becomes that of a perfect fluid with isotropic pressure. In this case po 
and Uq correspond to the pressure and energy density of a classical neutron gas. 



5. Local kinematic variables. 

Since we are interested on the local evolution of volume elements of the magnetized neutron 
gas associated with the source ( |3.4|) , we need to re-write the dynamical Einstein-Maxwell 
equations in terms of covariant parameters associated with the local kinematics of volume 
elements as described by u a . For the Kasner metric in the comoving geodesic 4-velocity 
frame, the nonzero local kinematic parameters are the expansion scalar, 0, and shear 
tensor, a a ^ , given by: 



6 = u% (5.1) 

e 
3" 



where h a p = u a up+g a p is the projection tensor and round brackets denote symmetrization 
on the indices a, f3. The geometric interpretation of these parameters is straightforward: 
denotes the isotropic proper time change of proper volume of local fluid elements, whereas 
a a a describes the deformation of local volumes as they expand at different rates along the 
directions given by its eigenvectors. 

The expansion scalar and components of the shear tensor for the Kasner metric are: 

VI W2 HZ 

a» v = diag[a%, a y y , a z z , 0] = diag [Ei, E 2 , S 3 , 0], (5.4) 



-9- 



where 

Sq = ItT ~ ^7T ~ ^7T' a^b^c, (a, b, c = 1,2,3). (5.5) 
3C^a 3C^f, 6LJ C 

Since the shear tensor is trace-free: a a a = 0, we can eliminate one of the quantities 
(Si, E2, S3) in terms of the other two. In fact, for the Bianchi I model in the Kasner 
metric, one of these quantities is enough to fully represent a a a , though for mathematical 
convenience we will keep two of these variables by eliminating Si in terms of (S2, S3). 
By means of ( |5.3j ) and (|5.5|), all second order derivatives of the metric functions in ([3.7|), 
(3.11) and ( 3.12j ) can be re-written as first order derivatives of 0, S2 and S3. After some 



algebraic manipulation we can re-write ( |3.7D , ( 3.11| ) and ( |3.12| ) as the following autonomous 
first order system of evolution equations 

U = -(U + p- -BM)9-BMZ 3 , (5.6) 
3 

■ kBM 

52 = BS 2 , (5.7) 

5 3 = 2 -kBM - 6S 3 , (5.8) 

6 = k{BM + ^{U -p))-Q 2 , (5.9) 

$ = ~/3(3E 3 -20). (5.10) 

Together with these equations, we have the following algebraic constraint 

B 2 

- S 2 - S 2 S 3 + — - S 2 = kU. (5.11) 



that follows from (3.10). This non-linear first order system in the variables U, /3, 0, £2, 



and the constraint (5.11) become fully determined once we use the thermodynamical equa- 
tions of the previous section to express M and B in terms of /3 = B/B c . The solution of 
this system describes the dynamical evolution of local volumes of a magnetized neutron 
gas that could be taken as a rough approximation to a grand canonical subsystem of this 
source near the center of a compact object. 

6. Dynamical equations. 



The system of evolution equations (5.6)-(p30|) can be transformed into a proper dynamical 



system by introducing the following dimension-less evolution parameter 

H = e ± = ±l, (6.1) 

3 dr Ho dt 
together with the dimension-less functions: 

* = !■, S 2 = ^, 5 3 = f , /3 = f, (6.2) 

Hq Hq Hq D c 

where Hq is a constant inverse length scale, which for convenience we choose as 3Hq = 
k\ =>• \Hq\ = 1.66x 10 -4 cm -1 . Note that Hq is not the cosmological Hubble constant, given 
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by Hq OSI11 = 0.59 x 10 _28 cm _1 , but a constant that provides a length scale 1/Hq ~ 6 km 
that is adequate for the characteristic length scale of the system under consideration. The 
functions £2 and S3 are the shear tensor component normalized to this scale, while the 
dimension-less time r can be positive or negative, depending on the sign of Hq 
Substituting the variables 



± v / k\/3. 

into the system ( |5.6| ) and the constraint ( |5.11| ) we get 



S2,r 



^— (2H-S 3 )(F M 

-pr M - 3S 2 n, 

2/3T M - 3S 3 H, 

3r P 1 



2/3(S 3 - ZH), 



S2S3 



1 



■XSl+S; 



3J) 



(6.3) 

(6.4) 
(6.5) 

(6.6) 
(6.7) 



- Sf - 5| - S 2 S 3 + 3H 2 = 3r v . (6.8) 

where T indicates derivative with respect to r. We have also replaced the variable U for 
fi, because from (|2|) we have: U = U(f3,fi) => U jT = A(rj/ )At // )T + Tjj^^ t ), which allows us 
to obtain /u jT from U T . We will solve numerically the system fl6,3|) in the following section. 

7. Numerical solutions and physical discussion 

Since we are interested in studying a collapsing magnetized neutron gas configuration, we 
need to consider the local collapse of volume elements. Hence, we solve the constraint ( |6.8| ) 
to obtain the two roots of H, so that the condition for a collapsing evolution follows by 
selecting the negative root. To ensure the local collapse we shall demand in the solution of 
equations ( |6.3| ) that the initial expansion, 0, is negative, which implies for an initial time 
hyersurface r = that H(0) < 0. This follows from (|5.3| ) and ( |6.2| ) expressed in terms of 
the local proper volume V = y/ det g a p = Q1Q2Q3 as: 



y = V(0)exp^3y Udrj. (7.1) 
To investigate the direction dependence of the collapse, in terms of the directions given by 



the space coordinates (x,y,z), we can relate by means of (|5.3|) , fl5.5| ), (|6.1| ) and (|6.2j ) the 
spatial metric components to the combination T-L + Sj, leading to: 

Qjir) = g,(0) exp[ J (H + Sj)dr], j = 1, 2, 3. (7.2) 
where Qj(0), are constants that can be identified with initial values for Q,(r). 



To solve the system (O) we shall use a wide range of initial conditions associated with 
typical conditions prevailing in a neutron star Jl], [2|, [3|, [|, ||, ||, 0, g|, for example: /j, = 2 
p ~ 10 15 g/cm 3 , /?o = — for magnetic fields between 10 18 gauss and 10 15 gauss. 

We shall impose in all numerical trials the condition of volume collapse: T-L(0) < 0, together 
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with: 52(0) = 0, ±1, #3(0) = 0, ±1, which correspond to cases with zero initial deformation 
and initial deformation in the direction of the axes y or z. 

The numeric solution for the function displayed in Figure [l] for different initial condi- 
tions, shows that H — > —00, regardless of the selected initial conditions. The magnetic 
field tends to increase, but always remains below the critical field B = B c . This behavior 
is shown in Figure [l] for the full range of initial conditions. The plots displayed in Figure 
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Figure 1: The left panel shows the behavior of H vs r for different initial conditions. The right panel 
shows the magnetic field intensity (/? = B/B c ), it has tendency to rise, but remains below the value of the 
critical field. 

|2] and in the left panel of Figure |3| describe the collapse of the fluid elements from the 
solutions of (|6.3|) . It is evident from these figures that the quantities Si + H. — > — oo, 
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Figure 2: Behavior of (Si + H) and (S2 + H) versus r. We can see this quantities tending to -co, besides 
different collapse times for different initial conditions. 

so that the spatial metric coefficients tend to zero (Qi,Q2,Q3 — ^ 0), which clearly shows 
that volume elements collapse to a point like isotropic singularity. 



8. Phase Space. 



As it was done in [ 10 1 , we can use the constraint ( |6.8| ) to transform the evolution system 
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Figure 3: The left panel shows behavior of (§3 + H) versus r. The behavior is similar to the trajectories 
given in Fig^, therefore the quantity (S3 + H) — > —00. The right panel shows paths in the small section of 
the space phase (S3,/3,/z). The points represent the unbranded initial conditions, the point "a" represents 
the attractor point which coordinates are (S3 = 0,/3 = 0, fj, = 1, H — 0). 



plf) to a reduced system of equations in the variables S3,/3,^, H 



1 



/': 



s 3 , T = 2/?r M - 3S 3 n 

H,r • 3 

At 



(ZH - S 3 )(T M - 2T Ut p)P - m{T P + Tu) 



2/3(5 3 - 2H), 



T P ) - m z 



(8.1) 
(8.2) 
(8.3) 
(8.4) 



where we note that the only equation that was modified is the equation for Ji, hence (| 
is equivalent to (|Q|)- 

The trajectories in the 3-dimensional subsection of the phase space, parametrized by 
(S3,{3,[i), are shown in Figure ||[ The evolution of the system is determined by the sign 
of H . For r < 

H 



Hq = — y^vX/3, the system evolves towards the stable attractor 
(point marked by a), while for r > 



yj AtA/3, the trajectories evolve towards a 
singularity. A similar study was conducted for the remaining 3-dimensional subsections of 
the space phase, obtaining qualitatively similar results, which was obtained coordinates of 
attractor are: (S3 = 0, f3 = 0, fi = 1, % = 0). 



9. Conclusions. 

We have used a Bianchi I model to study the evolution of a magnetized neutron gas 
characterized by a physically motivated and fully relativistic equation of state. As far 
as we are aware, this equation of state has not been considered previously in a general 
relativistic context. Besides the general theoretical interest in undertaking such a study, 
we argue that the simplified Bianchi geometry roughly approximates a grand canonical 
subsystem of a magnetized neutron source in the conditions prevailing near the center 
of a compact object, hence our treatment can be conceived as a toy model that can be 
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useful in understanding the local evolution of volume elements of this source under these 
conditions. However, a proper examination of the specific limitations of the dynamics of 
this toy model and/or its connection with concrete astrophysical studies of actual compact 
objects lies beyond the scope of the present article. As we comment further ahead, we 
will consider these important tasks in forthcoming articles by resorting to perturbation 
techniques, more elaborate numerical methods and less idealized sources. 

The Einstein-Maxwell field equations for a magnetized neutron gas in the Bianchi I 
geometry were transformed into a set of non-linear evolution equations, which were solved 
numerically for generic collapsing initial conditions and analyzed qualitatively as a proper 
dynamical system. The results that we found are: 

• The final state in the collapsing evolution of local volume elements is an isotropic 
point -like singularity. This final state occurred for a wide range of initial conditions 
associated with parameter values that would be typical in a compact object. 

• The magnetic field increases, but its values are always below the critical field. This 
result is consistent with numerical values of maximal field intensities compatible with 
stability in numerical studies of magnetized rotating configurations (see Q). 

• The study of the phase space associated with the dynamical equations shows that 
the system evolves, for Hq < 0, to an equilibrium point, (i.e. into a stable configura- 
tion). It is possible to introduce a temperature dependence in the equation of state. 
Buy doing so, the of evolution of the neutron gas could be associated with high tem- 
perature neutron sources in the context of early universe conditions in cosmological 
models dealing with primordial magnetic field |2^| . 

It is important to remark that, unlike the dynamical study of a magnetized electron 
gas [jl0|, anisotropic "cigar-like" singularities did not occur for a wide range of initial 
conditions. Since a dynamical effect of the magnetic field under critical conditions is a final 
state anisotropic singularity aligned in the direction of the field, the exclusive occurrence 
of point-like isotropic in a magnetized neutron gas suggests that the final stage of the 
evolution of this gas is more intensely dominated by the focusing effect associated with 
strong gravity than the electron gas. This is consistent with the fact that electrons are 
strongly coupled to the magnetic field through their electric charge, whereas neutrons have 
a weaker coupling associated with their anomalous magnetic moment. 

It is important to stress that the geometry of Bianchi I models in a comoving frame 
has stringent limitations in dealing adequately with the dynamical effects associated with 
a magnetic field. This is important when considering neutrons as a source in which electric 
charge vanishes but not the magnetic moment. By being spatially homogeneous with a 4- 
velocity orthogonal to flat 3-dimensional hypersurfaces of maximal symmetry, the Lorentz 
force is necessarily zero: f b = qu a F ab = 0. Also, the fact that the Bianchi I model 
is spatially flat makes it inadequate to examine (even as a toy model) the interplay of 
local collapse and the magnetic field in localized objects, as such interplay is necessarily 
associated with strong positive spatial curvature. However, in our case these inadequacy 
can be overcome (at least partially) by considering a general perturbation scheme on a 
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Bianchi I model, in which spatial curvature and 4-acceleration are perturbative but not 
zero. In this case it is possible to examine the effects of the magneto-curvature coupling 
associated with a non-trivial Lorentz force and a nonzero deceleration parameter in the 
Raychaudhuri equation (see |]l^] for general detail). The use of such a perturbed Bianchi 
I model for the description of the neutron gas source considered in this article is presently 
under consideration in a separate article. 

Besides the introduction of a perturbation scheme in a Bianchi I model, another possi- 
ble improvement on the dynamical description of the source under consideration would be 
to consider Bianchi models I, V, VII or IX with a tilted 4-velocity, which are endowed with 
more degrees of dynamical freedom, including even the possibility of nonzero rotation (see 
|23tl and references quoted therein) . These models would allow for a less restrictive study of 
the dynamical effects, reported in |i~8|| , in which magnetic tension and gravitational collapse 
may present non-trivial coupling with a nonzero and non-perturbative 4-acceleration and 
vorticity with the magnetic field. 

Finally, as we mentioned in the introduction, a magnetized gas consisting only of 
neutrons can be theoretically interesting but it is too idealized as a potential source for 
a compact object. Thus, we will consider as an extension of this work, besides the extra 
degrees of freedom in the dynamics (mentioned above), a gas mixture of neutrons, electrons 
and protons, complying with suitable balance conditions and adequate chemical potentials, 
in comparison with other types of equations of state (2|, ^, g]. These extension of the present 
work are also under consideration for future articles. 
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